Reciprocal transmittances and reflectances: An elementary proof 
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We present an elementary proof concerning reciprocal transmittances and reflectances. The proof 
is direct, simple, and valid for the diverse objects that can be absorptive and induce diffraction and 
scattering, as long as the objects respond linearly and locally to electromagnetic waves. The proof 
enables students who understand the basics of classical electromagnetics to grasp the physical basis 
of reciprocal optical responses. In addition, we show an example to demonstrate reciprocal response 
numerically and experimentally. 



I. INTRODUCTION 

Reciprocity, which was first found by Lorentz at the 
end of 19th century, has a long history^ and has been 
derived in several formalisms. There are two typical re- 
ciprocal configurations in optical responses as shown in 
Fig. [TJ The configurations in Figs. [Ha) and QJb) are 
transmission reciprocal and those in Figs. [Ha) and [He) 
are reflection reciprocal. As shown in Fig. [1] we denote 
transmittance by T and reflectance by R; the suffice k 
and stand for incident wavenumber vector and angle, 
respectively. The reciprocal configurations are obtained 
by symmetry operations on the incident light of the 
wavenumber vector: (k x , k z ) — > {—k x ,—k z ) or (— k x , k z ). 
Reciprocity on transmission means that 7k = T_k, and 
that on reflection is expressed as Re — R-e, which is not 
intuitively obvious and is frequently surprising to stu- 
dents. 

The most general proof was published by Petit in 
1980, 2 where reciprocal reflection as shown in Fig. Q] is 
derived for asymmetric gratings such as an echelette grat- 
ing. On the basis of the reciprocal relation for the solu- 
tions of the Helmholtz equation, the proof showed that 
reciprocal reflection holds for periodic objects irrespec- 
tive of absorption. It seems difficult to apply the proof to 
transmission because it would be necessary to construct 
solutions of Maxwell equations that satisfy the bound- 
ary conditions at the interfaces of the incident, grating, 
and transmitted layers. The history of the literature on 
reciprocal optical responses has been reviewed in Ref Q]. 

Since the 1950s, scattering problems regarding light, 
elementary particles, and so on have been addressed by 
using scattering matrix (S-matrix). In the studies em- 
ploying the S-matrix, it is assumed that there is no ab- 
sorption by the object. The assumption leads to the uni- 
tarity of the S-matrix and makes it possible to prove reci- 
procity. The reciprocal reflection of lossless objects was 
verified in this formalism. 3 

In this paper we present a simple, direct, and general 
derivation of the reciprocal optical responses for trans- 
mission and reflection relying only on classical electrody- 
namics. We start from the reciprocal theorem described 



in Sec. |TT]and derive the equation for zeroth order trans- 
mission and reflection coefficients in Sec. Mil The equa- 
tion is essential to the reciprocity. A numerical and ex- 
perimental example of reciprocity is presented in Sec. lIV 
The limitation and break down of reciprocal optical re- 
sponses are also discussed. 



II. RECIPROCAL THEOREM 

The reciprocal theorem has been proved in various 
fields, such as statistical mechanics, quantum mechanics, 
and electromagnetism. 4 Here we introduce the theorem 
for electromagnetism. 

When two currents exist as in Fig. [2] and the induced 
electromagnetic (EM) waves travel in linear and locally 
responding media in which Di(r) — ^2j£ijEj(r) and 

B i( r ) = Y., /'-,//, ii'K then 



J ji(r) • E 2 (r)rfr = J j 2 (r) • Ei(r)dr. 



(1) 



Equation |T]) is the reciprocal theorem in electromag- 
netism. The proof shown in Ref. |4| exploits plane waves 
and is straightforward. Equation ([1]) is valid even for 
media with losses. The integrands take non-zero values 
at the position r where currents exist, that is, jj(r) ^ 0. 
The theorem indicates the reciprocity between the two 
current sources ji (i — 1,2) and the induced EM waves 
Ej which are observed at the position of the other source 

j fe (M*)- 



III. RECIPROCAL OPTICAL RESPONSES 

In this section, we apply the reciprocal theorem to op- 
tical responses in both transmission and reflection con- 
figurations. First, we define the notation used in the 
calculations of the integrals in Eq. fl} . An electric dipolc 
oscillating at the frequency u> emits dipole radiation, 
which is detected in the far field. When a small dipolc 
p along the z axis is located at the origin, it is written 
as p(t) = p(t)e z and p(t) = poe tuJt , where e z denotes the 
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unit vector along the z axis and po the magnitude of the 
dipole. The dipole in vacuum emits radiation, which in 
the far field is 



E(r,t)= / ^ sin< 
v ' Atteq c 2 r 



eo 



-e smw-ee, 



47reo err 



(2a) 
(2b) 



where polar coordinates (r, 9, 4>) are used, a unit vector 
is given by e s = (cos 9 cos 0, cos sin 0, — sin 9) , and t' — 
t — r/c. Because the dipole p is defined by p(r, t) — 
j rp(r, t)dr and conservation of charge density is given 
by V • j + dp/dt — 0, we obtain the current j associated 
with the dipole p: 



j(r,f) = p(t)d(r)e z 



(3) 



Consider two arrays of N dipoles (long but finite) in 
the xz plane as shown in Fig. [3J The two arrays have the 
same length, and the directions are specified by normal- 
ized vectors (i — 1,2) and ni || 112. In this case, the 
current is ji || n^. If the dipoles coherently oscillate with 
the same phase, then the emitted electric fields are su- 
perimposed and form a wave front at a position far from 
the array in the xz plane as drawn in Fig. [3J The electric 
field vector of the wave front, E^j n , satisfies Ej^n || and 
travels with wavenumber vector k^. Thus, if we place 
the dipole arrays far enough from the object, the induced 
EM waves become slowly decaying incident plane waves 
in the xz plane to a good approximation. The arrays of 
dipoles have to be long enough to form the plane wave. 

For the transmission configuration, we calculate J ji • 
Efc dr (i,k =1,2 and i ^ k). Figure [3] shows a typical 
transmission configuration, which includes an arbitrary 
periodic object asymmetric along the z axis. The relation 
between the current , the direction of the dipole, and 
the wavenumber vector k^in of the wave front is summa- 
rized as ji || rii and _L ki ; ; n . It is convenient to expand 
the electric field into a Fourier series for the calculation 
of periodic sources: 



E(r)=^E( m )exp(zk m -r), 



(4) 



where E^" 1 ) is the Fourier coefficient of E(r), k m = 
[kx.m 1 0, k zm ) = [kin,x "F 2Trm/d x , 0, k zm ) (m = 
0, ±1, ±2, • • • ), and d x is the periodicity of the object 
along the x axis (see Fig. [3J. The z component is ex- 
pressed in homogeneous media in vacuum as k Zjm = 

±^Jlif n — k% rnl where the signs correspond to the direc- 
tions along the z axis. 

When the dipole array is composed of sufficiently small 
and numerous dipoles, the integration can be calculated 



to good accuracy as 

/ ji(r) • E 2 (r)dr = / icupo^ ■ ^ E^" l) exp(ik m • sn^ds 
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(5a) 
(5b) 

(5c) 



where E. 



(0) ^ , E (o) 



To ensure that the integration is 



proportional to 5 m> o, the array of dipoles has to be longer 
than L: 



L = (length of dipole) • q, 



(6) 



where q is the least common multiple of the diffraction 
channels which are open at the frequency uj. This condi- 
tion would usually be satisfied when Ei ; n forms a plane 
wave. 

By permutating 1 and 2 in Eq. (|5c(l . we obtain Jj2 • 

Eidr = iu}NpoE[°\ Equation (|5c| and the reciprocal 
theorem in Eq. Jl} lead to the equation 



E (0) = E (o) 

? (0) 



(7) 



Each electric vector E\ (i = 1, 2) is observed at the 
position r where there is another current jfe(r) (k ^ i). 
The integral in Eq. fTJ) is reduced to Eq. ([5c| which is 
expressed only by the zeroth components of the trans- 
mitted electric field. The reciprocity is thus independent 
of higher order harmonics, which are responsible for the 
modulated EM fields in structured objects. When there 
is no periodic object in Fig.[3j a similar relation holds: 
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The transmittance Ti is given by 



T; 



E 



(0) 



E 



no,(o) 



(8) 



(9) 



From Eqs. ([T])-©, we finally reach the reciprocal relation 
Ti=T 2 . 

The feature of the proof that Ti = T2 is independent 
of the detailed evaluation of and therefore makes 
the proof simple and general. The proof can be ex- 
tended to two-dimensional periodic structure by replac- 
ing the one-dimensional periodic structure in Fig. [3] by 
two-dimensional one. Although we have considered pe- 
riodic objects, the proof can also be extended to non- 
periodic objects. To do this extension, Eq. has to be 
expressed in the general form E(r) = J E(k) exp(ik-r)dk, 
and a more detailed calculation for J jj • E^dr is required. 
Reciprocity for transmission thus holds irrespective of ab- 
sorption, diffraction, and scattering by objects. 

In Fig. [3J the induced electric fields E; are polarized in 
the xz plane. The polarization is called TM polarization 
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in the terminology of waveguide theory and is also often 
called p polarization. For TE polarization (which is of- 
ten called s polarization) for which Ej has a polarization 
parallel to the y axis, the proof is similar to what we have 
described except that the dipoles are aligned along the y 
axis. 

Reciprocal reflection is also shown in a similar way. 
The configuration is depicted in Fig. 2J The two sources 
have to be located to satisfy the mirror symmetry about 
the z axis. The calculation of fit ■ E^dr leads to the 
reciprocal relation for reflectance Ri = R 2 . Note that 
£*o,(0) j n jjgjj nag ^ Q k e evama t ec j ^ replacing the 

periodic object by a perfect mirror. 

IV. NUMERICAL AND EXPERIMENTAL 
CONFIRMATION 

An example of reciprocal optical response is shown 
here. Figure 03a) displays the structure of the sample 
and reciprocal transmission configuration. The sample 
consists of periodic grooves etched in metallic films of 
Au and Cr on a quartz substrate. The periodicity is 
1200 nm, as indicated by the dotted lines in Fig. EJa). 
The unit cell has the structure of Au:air :Au:air = 3:1:4:5. 
The thickness of Au, Cr, and quartz is 40 nm, 5nm, and 
1 mm, respectively. The structure is obviously asymmet- 
ric about the z axis. The profile was modeled from an 
AFM image of the fabricated sample. 

Figure 0(b) shows our numerical results. The incident 
light has 9 = 10° and TM polarization (the electric vec- 
tor is in the xz plane). The numerical calculation was 
done with an improved S-matrix method^' - The permit- 
tivities of gold and chromium were taken from Refs. H 
and 0; the permittivity of quartz is well known to be 
2.13. In the numerical calculation, the incident light is 
taken to be a plane wave, and harmonics up to n = ±75 
in Eq. ([5]) were used, which is enough to obtain accurate 
optical responses. The result indicates that transmis- 
sion spectra (lower solid line) are numerically the same 
in the reciprocal configurations, while reflection (upper 
solid line) and absorption (dotted line) spectra show a 
definite difference. The absorption is plotted along the 
left axis. The difference implies that surface excitations 
are different on each side and absorb different numbers of 
photons. Nonetheless, the transmission spectra are the 
same for incident wavenumber vectors and k2 i i n - 

Experimental transmission spectra are shown in 
Fig. [5jc) and are consistent within experimental error. 
Reciprocity is thus confirmed both numerically and ex- 
perimentally. There have a few experiments on reciprocal 
transmission (see references in Ref. [l]). In comparison 



with these results, Fig. [5](c) shows the excellent agree- 
ment of reciprocal transmission and is the best available 
experimental evidence supporting reciprocity. 

We note that transmission spectra in Figs. [5fb) and 
M,c) agree quantitatively above 700 nm. On the other 
hand, they show a qualitative discrepancy below 700 nm. 
The result could come from the difference between the 
modeled profile in Fig. [5Ja) and the actual profile of the 
sample. The dip at 660 nm stems from a surface plasmon 
at the metal-air interface, so that the measured transmis- 
sion spectra would be affected significantly by the surface 
roughness and the deviation from the modeled structure. 



V. REMARKS AND SUMMARY 

As described in Sec. lU the reciprocal theorem assumes 
that all media are linear and show local response. Logi- 
cally, it can happen that the reciprocal optical responses 
do not hold for nonlinear or nonlocally responding media. 

Reference [13 discusses an explicit difference of the 
transmittance for a reciprocal configuration in a non- 
linear optical crystal of KNbC^Mn. The values of 
the transmittance deviate by a few tens of percent 
in the reciprocal configuration. The crystal has a 
second-order response such that Di(r) = J^j £ ijEj( r ) + 
k £ ijkEj(r)Ek(r). The break down of reciprocity 
comes from the nonlinearity. 

Does reciprocity also break down in nonlocal media? 
In nonlocal media the induction D is given by D(r) = 
Je(r, r')E(r')dr'. Although a general proof for this case 
has not been reported to our knowledge, it has been 
shown that reciprocity holds in a particular stratified 
structure composed of nonlocal medial 

In summary, we have presented an elementary and 
heuristic proof of the reciprocal optical responses for 
transmittance and reflectance. When the reciprocal the- 
orem in Eq. |l| holds, the reciprocal relations come from 
geometrical configurations of light sources and observa- 
tion points, and are independent of the details of the ob- 
jects. Transmission reciprocity has been confirmed both 
numerically and experimentally. 
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FIG. 1: Reciprocal configurations, (a) and (b) show recipro- 
cal configurations for transmission. Tk in (a) denotes trans- 
mittance for incident wavenumber vector k. T_k in (b) is 
defined similarly. The reciprocal relation is Tk = T_k- (a) 
and (c) are reciprocal for reflection. Re in (a) is reflectance 
for incident wavenumber vector (k x ,k z ) and R-e in (c) for 
(— k x ,k z ). The reciprocal relation is Rg — R~o- 
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FIG. 2: Schematic drawing of two currents ji and the electric 
fields E, induced by ji (i — 1,2). The curves denote the 
position where the currents exist. 



5 




FIG. 3: Schematic drawing of reciprocal configuration for 
transmission. The object has an arbitrary periodic structure, 
which is asymmetric along the z axis. Currents ji induce 
electric fields Ej,in (i = 1, 2). 




FIG. 4: Schematic configuration for reciprocal reflection. The 
object has an arbitrary periodic structure, which consists of 
asymmetric unit cells. The currents ji yield electric fields 
Ei, in (i = 1,2). 
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FIG. 5: (a) Schematic drawing of metallic grating profile mod- 
eled from AFM images. The periodicity is 1200 nm. The dot- 
ted lines show the unit cells in which the ratio is Au:air:Au:air 
= 3:1:4:5. The thickness of Au, Cr, and the quartz substrate 
is 40 nm, 5 nm, and 1 mm, respectively, (b) Numerically cal- 
culated spectra for 10° incidence of ki,j n (upper panel) and 
k2.in (lower panel) of TM polarization. In both panels the 
reflectance (upper solid line) and absorption (dotted line) are 
plotted using the left axis, while the transmittance (lower 
solid line) uses the right axis, (c) Measured transmittance 
spectra, corresponding to the transmittance spectra in (b). 



